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Abstract 

Using a completely analytic procedure - based on a suitable exten- 
sion of a classical method - we discuss an approach to the Poincare- 
Mel'nikov theory, which can be conveniently applied also to the case 
of non-hyperbolic critical points, and even if the critical point is lo- 
cated at the infinity. In this paper, we concentrate our attention on 
the latter case, and precisely on problems described by Kepler-like 
potentials in one or two degrees of freedom, in the presence of general 
time-dependent perturbations. We show that the appearance of chaos 
(possibly including Arnol'd diffusion) can be proved quite easily and 
in a direct way, without resorting to singular coordinate transforma- 
tions, such as the McGehee or blowing-up transformations. Natural 
examples are provided by the classical Gylden problem, originally pro- 
posed in celestial mechanics, but also of interest in different fields, and 
by the general 3-body problem in classical mechanics. 
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In a previous paper^, we discussed a completely analytic procedure - 
based on a suitable extension of a classical method^ - to introduce the 
Poincare-Melnikov theory^"^ concerning the appearance of chaotic behaviour, 
which can be applied even to the case where the critical point is not hyper- 
bolic. The case of nonhyperbolic stationary points has been already consid- 
ered by several authors, although with quite different methods or in different 
contexts: we mention here Ref.s 6-11; some other references, more strictly 
related to our arguments, will be given in the following. In Ref. 1, we have 
also shown that our procedure may work even in some problems where the 
critical point is located at the infinity of the real line, as occurs in the case 
of the classical Sitnikov problem^^. 

In this paper, we want first of all to reconsider and refine this approach, 
concentrating our attention precisely on problems described by Kepler-like 
potentials, in the presence of time-dependent perturbations of very general 
form. An example is provided by the classical Gylden problem originally 
proposed in celestial mechanics, but also of interest in different fields (see 
Ref. 13). We shall then extend this procedure to problems with 2 degrees 
of freedom and in the presence of rotational symmetry, and show that the 
appearance of chaos (possibly including Arnol'd diffusion^^) can be proved 
quite easily and in a direct way, without resorting to singular coordinate 
transformations, such as the McGehee transformations-*^^. 

Natural examples are provided by the general planar 3-body problem in 
classical mechanics^^"-*^^. 

Let us remark that, although our procedure is quite general, we shall 
consider here - for the sake of definiteness and clarity - only the case of the 
Kepler potential V{r) — 1/r: from the following discussion it will appear 
completely clear that the method can be easily applied - with suitable minor 
adjustments - to different classes of problems, for instance to problems where 
the Kcplcr potential is replaced by another "long range" potential, as e.g. 
V 1/r^. 
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I. The Gylden-like problems 



We start considering a 2-degrees of freedom problem described by a 
Hamiltonian of the form 

H=lp'-- + eW{r,t) (1) 

where p = (pi,P2), i-e a standard (planar) Kepler problem plus a smooth 
perturbation eVF(r, t). This problem can be easily reduced to a 2-dimensional 
dynamical system for the two variables r, r, and for this reason we shall treat 
first this case, not only for better illustrating our procedure, but also in view 
of some direct and interesting applications, which include the classical Gylden 
problem^^. 

Consider the "parabolic" solution of the unperturbed (e = 0) problem 
(1), which plays here the role of the homoclinic solution corresponding to 
the critical point at the infinity, i.e. r = oo, r = 0: it satisfies the equations 

^ = ± (^ = -2 (2) 

where A; 7^ is the (constant) angular momentum, and the sign — (resp. +) 
holds for t < (resp. t > 0). Prom (2) one gets 



±t — — - — V2r — k'^ + const 6 — ±2 arctan h const (3) 

3 k 

Let us denote by 

R = R{t) and Q = Q{t) (4) 

the expressions giving the dependence of r and of 6 on the time t which are 
obtained "inverting" equations (3) with the conditions R{0) — Tmin — 
and 0(0) = TT (let us emphasize that it is not necessary, for our purposes, to 
have the explicit form of these functions). It will be useful only to remark 
that R{t) is an even function and G(t) is a odd function of time. The choice 
6(0) = TT corresponds to select the solution describing the parabola with 
axis coinciding with the Xi axis and going to infinity when xi — > +00; the 
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whole Hamiltonian (1), including the perturbation W{r,t), is rotationally 
invariant, and therefore this choice is clearly not restrictive. We refer to the 
next section, where we shall consider symmetry-breaking perturbations, for 
some comment on this aspect. 

The first point is now to find the negatively and positively asymptotic 
sets to the critical point r — oo, r = in the presence of the perturbation 
eW{r, t); notice that this point would in fact correspond to the critical point 
X — y — under the McGehee's singular coordinate transformation^^'^^'^^ 

r — — r = y dt — — ds (5) 

x-^ 

but - instead of using this transformation (which usually requires quite cum- 
bersome calculations), we introduce a direct method similar to the classical 
one used in Ref. 2 (see also Ref. 1). Precisely, we want to show that some 
natural assumptions on the perturbation W{r, t) may guarantee not only 
the existence of smooth solutions approaching r = cxd, r = for t ±cxd, 
playing in this context the role of stable and unstable manifolds of the crit- 
ical point at the infinity, but also the possible presence of infinitely many 
intersections of these asymptotic sets on the Poincarc sections. 

It is important to remark that the critical point at the infinity is clearly 
not a hyperbolic point, and therefore the standard results of perturbation 
theory valid for hyperbolic points cannot be applied here. In particular, in 
order to be granted that the perturbation will preserve the property of the 
point at the infinity of being a critical equilibrium, we have to impose that 
the perturbation vanishes at the infinity. We shall give the precise rate of 
this vanishing in the following (see eq. (18) below). 

Let us write problem (1) in the form 

and let us look for solutions r{t) of the perturbed problem (6) "near" the 
family of the homochnic orbits R{t — to): we then put (see Ref. 2) 

r = R{t- to) + z{t - to) (7) 
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Inserting into (6) (with the time shift t — to ^ t), we obtain 

^+{^-^)' = ^^'^^^^^ + ^o^ (S) 

where the r.h.s. (which we will shortly denote by G{t, to)) is given by 

dW(R(t).t + to) , , . / X 

G{t, to) = — e ^ — — + higher order terms m z[t) (9) 

or 

Consider the homogeneous equation obtained putting e = in (8): one 
solution is clearly R{t), another independent solution ilj{t) can be constructed 
with standard methods (see e.g. Ref. 19) or by direct substitution: these 
two solutions have a different behaviour for t — > ±oo, precisely: 

R(t) ~ \t\-^/^ and ijjit) ~ \t\^/^ oo (10) 

As well known, the general solution z{t) of the complete (nonhomogcncous) 
equation (8) can be written in the following integral form, with A, B arbitrary 
constants and ti arbitrarily fixed^'^^: 

z{t) ^ AR{t)+Bil){t)-R{t) ilj{s)G{s,to) ds+ilj{t) R{s)G{s,to) ds (11) 

Jti Jti 

Let us now look for solutions z^~\t) (and resp. z^~^\t)) of (11) with the 
property of being bounded for t — > — oo (resp. t — > +oo); these will provide 
solutions 

rW(i) = i? + ^W 

of (6) which belong, by definition, to the unstable and stable manifolds of 
the point r = oo, r — 0. From (11), we have then to require that the two 
quantities 

i){t)(^B+J^'R{s)G{s,to)ds) and R{t)(^A- J^' '4){s)G{s,to) ds) (12) 

remain bounded as t — > — oo when looking for the r^~^ solutions, and resp. 
as t — > +00 for the r^+^ solutions. 
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R(t)l'^(sf-Wll±lAi, (14) 



Consider now the linearization of the problem (11) around the solution 
z{t) = 0: this amounts in particular to deleting the higher-order terms in 
the expression of G in (9). Taking also into account the different behaviour 
given in (10) of the two fundamental solutions R{t) and ■ip{t), it is easy to see 
that the above conditions on the quantities (12) are simultaneously satisfied 
both Bit t — — oo and a,t t — +oo if for some to the following Mel'nikov-type 
condition 

p M).i + t„) ^ (13) 
j-oo or 
is fulfilled, together with the additional one, i.e. that the quantity 

dW{R{s),s + to) 

'ti 

is bounded as t ^ ±oo. Once these conditions are satisfied, one can conclude 
that, as a consequence of the implicit- function theorem^ (or also thanks to a 
suitable version of the Lyapunov- Schmidt procedure, see e.g. Ref. 20), there 
exists a smooth and bounded solution of (8). 

Let us now discuss the two above conditions (13-14). 

The first one is identical to the usual Mel'nikov condition obtained under 
the standard assumption that the critical point is hyperbolic^"^: eq.(13) can 
therefore be viewed as an extension of the classical Mel'nikov formula to the 
present "degenerate" case, in which the critical point is at the infinity. Let 
us now assume that the perturbation W{r,t) is a smooth function, periodic 
in the time t, with arbitrary period T (it is not restrictive to assume T — 27r) 
and zero mean-valued: 
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27r 

W{r,t)dt^O (15) 





Then one can consider its Fourier expansion 

oo 

W{r, t)^Yl i^n{r) cos nt + Bn{r) sin nt) (16) 

n=l 

and (thanks to the parity of the function R{t)) write down the Mel'nikov 
function M{to) defined in (13) in the form: 

oo 

-^(^o) = X] ("^^ "-^0 + Pn sin nto) (17) 
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where 

f+'^ ,dAJR(t)) . , ^ f+°° ,dBJR(t)) . 

R(t) — \ ^ sinnt dt , = / R{t) — , sinnt dt 
J —oo dv J —oo dv 

(17') 

Prom this expression, one can immediately conclude that the function M(io), 
being a smooth periodic function with vanishing mean value, must certainly 
possess zeroes, thus fulfilling condition (13). 

For what concerns the second condition, which requires the boundedness 
of (14) and which appears here to compensate the lack of the "exponential di- 
chotomy" peculiar of the hyperbolic case, a simple estimate of the behaviour 
as |t| — >■ cx) of the integral in (14), using (10) and (16), and recalling from (3) 
that R{t) ~ |tp/^ as |t| oo, shows that it is sufficient to assume that the 
quantities An{r), Bn{r) in the expansion (16) vanish as r — > oo according to 

A(r)-^ , i?n(r)~^ with 5>l/2 (18) 

in order to be granted that the above condition on (14) is satisfied. It can 
be noticed that the same condition (18) would also guarantee that under the 
McGehee transformation^^'^^'^^ (5), the perturbation is not singular at x = 0. 

Changing now the point of view, and considering the Poincarc sections of 
the r^~^ and r*^+^ solutions, the above arguments show that, once conditions 
(13-14) are satisfied, there occurs a crossing of the negatively and positively 
asymptotic sets on the Poincare section. One usually imposes at this point 
that the zeroes of the Mel'nikov function (13) are simple zeroes, i.e. that 

which ensures the transversality of the intersections, recalling that the func- 
tion M(to) expresses the signed distance between the intersecting manifolds^'^. 
Actually, it can be noted that it is not strictly necessary to impose this con- 
dition, indeed - according to an interesting and useful result by Burns and 
Weiss^^ - it is sufficient that the crossing is "topological", i.e., roughly, that 
there is really a "crossing" , from one side to the other. But in our case this is 
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certainly satisfied, because M(to), being a smooth periodic and zero mean- 
valued function, must necessarily change sign (see also Ref. 22 for a careful 
discussion on non-transversal crossings). Using then standard arguments, 
which are not based on hyperbolicity, thanks to the periodicity of the per- 
turbation, one immediately deduces^'^'^'^^ that there is an infinite sequence of 
intersections, leading to a situation similar to the usual chain of homoclinic 
intersections typical of the homoclinic chaos. 

The presence of such infinitely many intersections is clearly reminiscent 
of the chaotic behaviour expressed by the Birkhoff-Smale theorem in terms of 
the equivalence to a symbolic dynamics described by the Smale horseshoes. 
Actually, this theorem cannot be directly used in the present context be- 
cause its standard proof is intrinsically based on hyperbolicity properties'^"''. 
However, several arguments can be invoked even in the present "degenerate" 
situation. First of all, for the case of degenerate critical points at the infinity, 
we can refer to the classical arguments used in Ref. 12, and reconsidered by 
many others (see e.g. Ref.s 13,16,23). More specifically, see Ref. 24, where 
an equivalence to a "nonhyperbolic horseshoe" has been proved, in which the 
contracting and expanding actions are not exponential but "polynomial" in 
time. Let us also notice, incidentally, that the presence of Smale horseshoes 
and of a positive topological entropy has been proved by means of a quite 
general geometrical or "topological" procedure^^ which holds, in the presence 
of area-preserving perturbations, even in cases of nonhyperbolic equilibrium 
points (i.e. not only in the case of degenerate critical points at the infinity, 
see Ref.s 1,25). Alternatively, in the general situation, one may possibly re- 
sort to the method of "blowing-up", devised to investigate the properties of 
nonhyperbolic singularities by means of suitable changes of coordinates^^'^''. 

Finally, for what concerns the regularity of the solutions and of the asymp- 
totic sets, see in particular Ref.s 23-25,28. 

Summarizing, we can state the following: 

Proposition 1. Consider a Kepler-like problem as in (1), and assume that 
the perturbation W{r, t) is a smooth time-periodic function with zero mean 
value (16). Assume that it vanishes with r — > oo m such a way that (18) are 
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satisfied. Then there is a chaotic behaviour of the solution, induced by a chain 
of infinitely many intersections in the Poincare section of the negatively and 
positively asymptotic sets of the critical point at the infinity. 

Let us conclude this section with the obvious remark that the case of the 
Gylden problem, for which the perturbation is given by 

W{r,t)^^ (20) 

where ii{t) is a periodic function, satisfies all the above assumption and 
therefore exhibits chaotic behaviour^^. The above discussion then generahzes 
this result to a larger class of problems and under weaker assumptions. 



II. Problems with 2 degrees of freedom. 

We now consider the case of planar Kepler-like problems as in (1) but in 
the presence of perturbations of the more general form W — W{r,9,t) 

= ^p' - - + eW{r, e, t) = Ho + eW{r, 9, t) (21) 
2 r 

In this case the reduction of the problem as performed in sect. I is no longer 
possible, and we have to handle the four variables xi,pi, X2,P2 (or r, r, 9, 9). 
The first point to be remarked is that the degeneracy of the critical point 
at the infinity, r = oo, r = 0, appears now even worse than before, indeed 
we have here a "continuous family of points at the infinity" , due to the ar- 
bitrarity of the angle 9; more precisely, the homoclinic manifold, i.e. the 
set of solutions of the unperturbed equation which are doubly asymptotic to 
r = oo, r = 0, is given here, for each fixed value A; 7^ of the angular mo- 
mentum, by the 2-dimensional manifold described by the family of parabolas 
of equations R = R{t — to), 0(t — to) + ^O) where R(t), 0(t) have been de- 
fined in sect. I (see (3-4)), with arbitrary to, 6*0, or - in cartesian coordinates 
u = {xi,pi,X2,P2) ~ by 

X = x{9o,t-to)= (22) 
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[R{t -to) COS {&{t- to) + eo), 

R(t - to) cos (0(t - to) + ^o) - R{t - to)Q{t - to) sin {Q{t - to) + ^o), 
R{t - to) sin {e{t- to) + 60), 

R{t - to) sin {e{t - to) + 9o) + R{t - to)e{t - to) cos {e{t - to) + 9o)) 

In order to find conditions ensuring the occurrence of intersections of stable 
and unstable manifolds for the perturbed problem, we follow a similar (suit- 
ably extended) procedure as in sect. I. We first look for smooth solutions 
near the homoclinic manifold (see sect. I; here, clearly, z = {zi, Z2, Z3, Z4)) 

u = x{Oo,t-to) + z{eo,t-to) (23) 

of the problem, which we now write in the form 

ii = JVuH = F{u) + eJVuW (24) 

(J being the standard symplectic matrix). Linearizing the problem along an 
arbitrarily fixed solution x(^o, ^ — ^o) in the family (22), we get the following 
equation for z{t) 

z^A{0o,t)z + G{0o,t,to) (25) 

where 

A{eo,t) = {VuF){xiOo,t)) (26) 

and 

^(^0, t, to) = eJ{VuW){R{t), e(t) + 9o,t + to) (27) 

All the solutions z(t) of (25) are are given by (cf. Ref. 19; when not essential, 
the dependence on 6*0, to will be sometimes dropped, for notational simplicity) 

z{t) = Zh{t) + $(t) /* <^-\s)G{s)ds (28) 
Jti 

where Zh{t) is any solution of the homogeneous linear problem 

z^A{9o,t)z (29) 
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and $ is a fundamental matrix of solutions of (29). There are certainly two 
solutions of (29), which are bounded for any t e R (and vanish for t — > icxo), 
namely 

I - % (30) 

as one may easily verify (this also follows from general arguments^^'^°). As 
seen in sect. I, due to the degeneracy of the critical point, instead of the 
exponential dichotomy, typical of the standard hyperbolic case, we now get a 
power behaviour |t|°" of the solutions, but the general arguments for control- 
ling the behaviour for large |t| of the solutions z{t) in (28) can still be used. 
Precisely (cf. Ref. 29), observing also that, for Hamiltonian problems, the 
matrix $~^J is the transposed of a fundamental matrix of solutions of the 
same problem (29), one deduces from (28) and (30) that there exist bounded 
solutions of (25), both for t — > +00 and for t — > —00, if the two following 
conditions are satisfied (cf. Ref.s 29,30) 



+°°^ ax(6'o,t) ^v^|y(^(^)^e(t) + eo,i + io))dt = Mi(^o,^o) = (31) 



dt 

and 



+00 i^MM , VuW{R{t), e{t) + 9o,t + to))dt = M2{9o, to) = (32) 
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where ( " ? " ) stands for the scalar product in R^. Proceeding just as in 
sect. I, we assume that the perturbation W{r, 9, t) is a smooth time-periodic 
function, and we still assume to hold a condition analogous to (18) in order to 
to guarantee the boundcdness of z{t) at t = ±00, as already discussed. The 
above conditions (31-32) can be more conveniently rewritten in the following 
form, using (22), 

M,{9o,to)= (33) 



f 

J — ( 



j^^^^ dWjRjt), 9(t) + 9o,t + to) ^ ^^^^dWjRjt), e{t) + 9o,t + to) 



dr 89 
dW{R{t),e{t) + Oo,t + to) 



di = 



M.(.., = r ''"-W')."^ +^..' + '.) ^ (34) 

J —00 Ou 
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It can be significant to remark here that it is easy to verify that these condi- 
tions are identical to the Mel'nikov conditions for the appearance of homo- 
clinic intersections given e.g. in Ref.s 5,31 in the standard hyperbolic case 
and deduced by means of a different procedure, namely 

(VuHo,JVuWYx(Oo,t),t + to)dt^ / {Ho,W}(x(---))dt^O (35) 

-oo ^ ' J —oo 

and 

(V^K,JVuW){x(Oo,t),t + to) dt^ / {K,W}(x{...)) dt^O (36) 

-oo J —oo 

where Hq is the unperturbed Hamiltonian, and K = X2P1—X1P2 is the angular 
momentum (which is indeed a constant of the motion for Hq). We can then 
say that, again and apart from the additional condition (18), just as in sect. I, 
our procedure provides an extension of these formulas to the nonhyperbolic 
degenerate case we are considering here. 

Notice also that, as a consequence of the vanishing of the perturbation W 
at t — ±00, the first Mel'nikov condition (33) can be written in the simpler 
form 

J -00 dt 
where clearly the derivative d/dt must be performed only with respect to the 
explicit time-dependence of W . 

It is also clear that, when the perturbation is independent of 6', as in the 
cases considered in sect. I, one consistently gets that the first condition (33) 
becomes just (13), whereas the second one (34) is identically satisfied. 

Let us now introduce the "Mel'nikov potential" W = >V(^o,^o) (cf. Ref. 
32), corresponding to the perturbation W{r, 9, t): 

/+00 
W{R{t),e{t) + eo,t + to)dt (37) 
-00 

then one gets from this definition and from (33'-34) 

dW dW 
Mi(^o, ^o) = ^ = M^iOo, to) = ^^0 (38) 

Oto OVq 
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In other words, the two Mel'nikov conditions are equivalent to the existence 
of stationary points for the Mel'nikov potential VV(^^O)^o)- On the other 
hand, W is a smooth doubly-periodic function, and such a function certainly 
possesses points 9o,to where the two partial derivatives in (38) vanish, and 
this imphes that the two conditions (33'-34) are certainly satisfied. 

Now, exactly the same arguments (and with analogous remarks) presented 
in sect. I show that the vanishing of the Mel'nikov functions entails the 
presence of a complicated dynamics, produced by the chain of the infinitely 
many intersections of the asymptotic sets: see, e.g., Ref.s 5,28 for a detailed 
description of this "multidimensional" case. 

In this context, we can also consider the appearance of Arnol'd diffusion^^'^ 
Indeed, the integrals of the Poisson brackets appearing in the Mel'nikov con- 
dition (35-36) give precisely the total amount of the "variations" produced 
by the perturbation, from t = —oo to t = +oo, to the quantities Hq and K 
along the homoclinic solution x- Proceeding in a similar way as in Ref. 14, 
we can now look for the intersections of the asymptotic sets corresponding 
respectively for t — > — oo and for t — > -|-oo to different values ki and k2 of the 
angular momentum K. Observing that for both these solutions the energy 
Hq is the same, Hq — O, then - following a by now classical idea^^ - the above 
Mel'nikov conditions must be replaced by conditions of the form 

Mi{9o, to) = M^iOo, to) + ki-k2^0 (39) 

and, exactly as in Ref. 14, the conclusion is that, for \ki — k2\ small enough, 
intersections occur even from homoclinic solutions corresponding to different 
values of K. This argument is then consistent with the occurence of Arnol'd 
diffusion: actually, a complete argument would necessitate a consideration 
of the role played, in the present "degenerate" situation, by the usual "non- 
resonance" conditions: anyway, we point out that, in the same situation, i.e. 
in the case of general 3-body problem, but using a different approach, the 
occurrence of Arnol'd diffusion has been discussed in great detail by Xia^''^^^. 
On the other hand, it is also known that the phenomenon of Arnol'd diffusion 
requires a quite delicate treatment: see Ref. 34 for a careful and updated 
discussion on this point. 
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Let us remark incidentally that the introduction of the above Mel'nikov 
potential W would be in general impossible if the perturbation is not Hamil- 
tonian, i.e. if the problem (24) has the form of a general dynamical system 

u = JVuHo + eg(u, t) = F(u) + eg(u, t) (40) 

where the perturbing term g{u,t) (or g{r,r,6,6,t)) is "generic". Then in 
this case the above arguments cannot be repeated, and in particular the two 
Mel'nikov conditions, which can now be written in the general form 

Mi(^o,io)= {VuHo,g){x{Oo:t),t + to) dt^O (41) 

M2{eo,to) = (VuK,g){x{Oo,t),t + to) dt ^ (42) 

give two "unrelated" restrictions on Q^.t^^. and one then remains with the 
problem of discovering if there are or not some which satisfy simulta- 

neously both these conditions. 

We can then summarize the above discussion in the following form. 

Proposition 2. Let us consider a perturbed Kepler problem with Hamil- 
tonian (21) where W is a smooth, time-periodic function, vanishing at the 
infinity according to (18). Then, there is a chaotic behaviour (possibly giving 
rise also to Arnol'd diffusion), induced by an infinite sequence of intersections 
in the Poincare section of the negatively and positively asymptotic sets of the 
critical point at the infinity. In the case where the perturbed problem has 
the form (40) with a non-Hamiltonian perturbation g{u,t), the same result 
is true if there are some 9Q,to satisfying simultaneously the two conditions 

(41-42). 
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III. Some applications and final remarks 



We shortly consider here some applications of the discussion presented in 
sect. II. 

As a first particular case, assume that the perturbation W{r,6,t) is of 
the form 

W^W{r, aO + bt) (43) 

where a, b are arbitrary constants (the constant a should be clearly an inte- 
ger, and & 7^ 0), then the two conditions (33'-34) actually coincide; observing 
on the other hand that the Fourier expansion of such a is a series as in 
(16) in terms of the single variable a6 + bt, then these two conditions take 
the same form as in (17), and the Mel'nikov function is a smooth periodic 
function of aO^ + bto, leading thus directly to the same conclusions obtained 
above for what concerns the existence of zeroes, and of their properties as 
well (it is really not a restriction to assume that W is zero mean-valued, cf. 
Ref. 13). 

A specially important example of this situation is provided by the re- 
stricted circular 3-body problem, in this case indeed the perturbation is given 
byi6 

W^'--'^^'-'^- , ' (44) 

^l + r2 + 2rcos(e-t) 

and therefore just one condition has to be considered. The presence of the 
chaos produced by the chain of intersections of the asymptotic sets then is 
automatically granted by our discussion. Notice that the above expression 
(44) is actually the first-order expansion of the full potential in terms of the 
parameter e (which in this case is given by the mass ratio ^ between two 
celestial bodies), but also the exact expression of this potential, as given in 
Ref. 16, is in fact a function oi 9 — t only. 

It is completely clear that, in the presence of a more general perturbation, 
e.g. of the form (just to give an example) 

W{r, e, t) = 1^(1) (r, e-t)+ W^2) (r, 9 + t) 
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one now obtains two different Mel'nilcov conditions, and the presence of a 
cliaotic beliaviour follows from the existence of simultaneous solutions 60, to, 
which is ensured by the arguments shown in sect. II. 

The above results hold essentially unchanged if the perturbation depends 
on two (or more) parameters ei, 62, . . ., i.e. if one assumes that W may be 
written in the form (cf. Ref. 29) 

W = eiW^(i)(r, e, t) + e2VF(2)(r, 6,1) + ... (45) 

A natural example is provided by the planar 3-body problem, where one has 
to deal, in the more general elliptic case, with a quite complicated expres- 
sion of the perturbation containing three parameters (different masses and 
eccentricity^^"^^). Let us remark, however, that, at least in the simpler case 
of restricted elliptic problem, in which one has two parameters (ei — /lis one 
mass ratio and €2 = 6 the eccentricity), the perturbation cannot be written 
as in the above "first-order" form (45), but rather it takes the form^^ 

W = /.(W^(i) {r,9,t) + e W^2) (r, 9,t)) (46) 

in which the eccentricity plays the role of a "second-order" perturbation. 

On the basis of our previous arguments, we can just say that the presence 
of zeroes of the Mel'nikov functions, as obtained above for the circular case 
62 = e = 0, cannot be destroyed by the higher-order perturbation due to the 
eccentricity, and therefore chaos should be expected to persist in the elliptic 
case, whereas Arnol'd diffusion should appear as a second-order effect. A 
complete study of the general 3-body problem, and a full discussion of its 
chaotic properties, including the appearance of Arnol'd diffusion, together 
with several other dynamical features, is given in Ref.s 16-18. 
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